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The vacuum effective potential of three-dimensional SU(2) lattice gauge theory 

in an applied color- magnetic field is computed over a wide range of field strengths. 

The background field is induced by an external current, as in continuum field theory. 

pg ' Scaling and finite volume effects are analyzed systematically. The first evidence from 

lattice simulations is obtained of the existence of a nontrivial minimum in the effective 

potential. This supports a "ferromagnetic" picture of gluon condensation, proposed 

60 . by Savvidy on the basis of a one-loop calculation in (3-|-l)-dimensional QCD. 
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The vacuum structure of quantum chromodynamics is known to play a fundamental 
role in strong interaction dynamics [|ll , yet a clear understanding of the physics underlying 
such basic properties of the vacuum as the gluon condensate is lacking. An attractive 
physical picture of gluon condensation was proposed more than ten years ago by Savvidy, 
who calculated the one-loop effective potential of the QCD vacuum in response to an applied 
color-magnetic field 0. The absolute minimum of the effective potential was found to lie 
at a nonzero value of the applied field, suggesting that the QCD vacuum behaves like a 
ferromagnet, with a condensate in the color-field arising spontaneously. This result was 
further analyzed by Matinyan and Savvidy |^, and was obtained independently by Pagels 
and Tomboulis Q. 

However, the one-loop calculation is not adequate a priori to support the ferromagnetic 
picture of gluon condensation, since perturbation theory becomes untrustworthy near the 
minimum of the effective potential. Moreover it was discovered by Nielsen and Olesen 
that a constant color- magnetic field generates unstable modes at the one-loop level [Q. 
Nevertheless, further analysis by the Copenhagen group suggested that the instability is 
removed by a nonperturbative structure that leaves the nontrivial minimum of the one-loop 
effective potential intact. [This approach was challenged more recently in Ref. however, 
where it was argued that the external field problem is fully nonperturbative, even in the 
region of large fields.] From another line of argument, Adler reached the conclusion that the 
existence of a truly nonperturbative minimum in the effective potential can be established 
from the leading-order calculation ^ . An attractive model of hadronic structure including 
the automatic formation of flux-tubes and bag-like excitations has also been developed based 
on the occurrence of a nontrivial minimum in the effective potential 0,0. 

In order to truly substantiate the ferromagnetic picture of the QCD vacuum, a genuine 
nonperturbative calculation is required. In principle, the effective potential is well suited 
to analysis in lattice gauge theory. In fact, two groups have recently attempted such a 
lattice calculation in pure gauge SU(2) [|T^,|rT|]. Unfortunately, the qualitative features of the 
effective potential obtained by the two groups are in sharp disagreement. They acknowledge 
significant nonscaling and finite volume effects in their calculations. Moreover, they were 
unable to draw definitive conclusions about the behavior of the effective potential at small 
fields, the region of greatest interest. 

We also note that very different algorithms for the introduction of an external field on 
the lattice were used in Refs. [|l^ and ]ll|, and we think that there are significant drawbacks 



to both methods. In Ref. [|r^ a link variable is introduced to account for the external field; 
the external link is aperiodic (corresponding to the usual external vector potential in the 
continuum theory), while periodic boundary conditions for the "dynamical" links are used. 
We think that this procedure leads to undesired interaction terms in the action on the 
lattice boundary. In Ref. [ill] on the other hand an Abelian projection technique is used to 
construct the lattice action, and the results for the effective potential turn out to depend 
strongly on the choice of operator that is diagonalized. 

In this Letter we present a calculation of the effective potential in three-dimensional 
pure-gauge lattice SU(2). Our rationale for doing the lattice calculation in three dimensions 
is two-fold. First, one of us has recently demonstrated that the one-loop effective potential in 
2+1 dimensions has the same qualitative features as in 3+1 dimensions |T^. In particular. 



the effective potential has its absolute minimum at a nonzero value of the applied color- 



magnetic field in both theories. Thus the Savvidy "ferromagnetic vacuum" is found to occur 
at the one-loop level in both (2+1)- and (3-1-1 )-dimensional QCD (the limitations of the 
one-loop calculation are of course the same in both cases). Secondly, one can more readily 
perform a thorough check of scaling and finite volume effects in a three dimensional lattice 
calculation that can be achieved (with the same computing power) for the corresponding 
four-dimensional system. 

We have also developed an algorithm for the introduction of a background field in non- 
Abelian lattice gauge theory (in an arbitrary number of dimensions) that is free of the 
drawbacks in the methods used in Refs. |]TD|,rT|, outlined above. Our algorithm is a straight- 
forward transcription of the procedure that is used in continuum field theory to induce a 
background field by a coupling to an external current. This approach has previously been 
considered in the context of the lattice Abelian Higgs model by Damgaard and Heller |13 



In the continuum U(l) theory, the (Euclidean) action for a coupling of the gauge field to 
an external current is 
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where for our purposes the external current has the form 
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If F^^{x) is constant within a finite region, then jj^^*(a;) describes a solenoidal current which 
circulates around the boundary. An integration by parts in Eq. (^ leads to 
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where now the region of support of the external field is unrestricted. The path integral in 
this trivial theory leads to a "classical" behavior for expectation values such as {Ffj_^) = F' 



ext 



An extension of the above procedure to nonAbelian theories was proposed in Ref. pT 
The action is 
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where the current is a covariant generalization of Eq. 
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with T>'r{A) = d^S — gf A^ the usual covariant derivative. -F^^*'° is the nonAbelian 
field strength constructed out of A'f*''^. Following the one-loop calculations of the effective 
potential [pi^,|6|,p!2| , we consider an external field with a fixed direction in gauge space 
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An integration by parts in Eq. @) then yields 
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where 



^f;=3 ^ d,AT' - d^Air'- (8) 



Notice that the interaction term in Eq. (|^ only contains the derivative part ^F of the field 
strength for the "dynamical" gauge fields A"^ (compare ^F with the full field strength tensor 
F^^ = dfj^A'^ — dyA"^ — gf^-^^'A^^Al). If the full field strength appeared in the interaction 
term, then Eq. (|^ would describe the same "trivial" physics as the U(l) action of Eq. 
(0). This difference between the two actions expresses the fact that the "charge" carried by 
nonAbelian gauge fields drives the different physics. 

Equations (^ and (|^ represent a sensible way to introduce an external field in the 
nonAbelian theory since they lead, at the classical level, to an "induced" gauge field that is 
equal to the applied field; moreover, at the one-loop level, this action leads to the Savvidy 
effective potential of interest. 

We now consider the lattice transcriptions of the above continuum theories. We specialize 
to three dimensions. The lattice equivalent of the U(l) action Eq. (|^) is [^ 
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where Sw is Wilson's plaquette action, the external field in lattice units F^^* is given by 

F-* = a'/'F^:\ (9b) 

and the "dynamical" field strength F^^^, is given as usual in terms of the imaginary part of 
the plaquette f/^,^: 

F^y = JpimU^,, P=^ (9c) 

(a the lattice spacing). In the SU(2) lattice theory: 

'5su(2) =<Sw -lYl ^ilu^ ^/^^*' (10a) 

sites 

where Sw is Wilson's plaquette action for SU(2), and 

UV^'^ 1 4 

^r = (^) ^^;f. ^-^. (lOb) 

To extract the derivative part ^F of the "dynamical" field strength on the lattice, we explicitly 
compute the commutator of any two links U^, Uy which span the plaquette: 

^F;=3 = -zy^Tr [{U^, - [f/^, f/J) \a^] , (10c) 

where a^ is the third Pauli matrix. The connection between the right hand side of Eq. ( |10c|) 
and ^F in the continuum limit follows from the identification of the link variables in this 
limit as U^ = exp(ia(7A"cr'^/2); the above lattice action reduces to the continuum action Eq. 
(0), up to corrections of 0{a), which is the same order of accuracy as the continuum limit of 



the Wilson action. [The couphng to the external current which leads to Eq. ( p.Oa|) breaks the 
local SU(2) gauge invariance of the Wilson action to a local U(l) symmetry. This implies, 
for example, that Tr(t/^j,cr^) is an invariant quantity in the theory defined by Eq. ( |10a| ).] Our 
construction of the lattice action using Eq. (|10c|) differs from the procedure followed in Ref. 
where ^F was defined by an Abelian projection (the results for the effective potential 
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obtained in Ref. [^ depend strongly on the choice of operator that is diagonalized) . 

In our calculations we take the field strength tensor for the external magnetic field 
scalar in three dimensions) to have nonvanishing components only in the (1,2) plane: 
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An important aspect of our lattice simulation is our use of free boundary conditions for 
the "dynamical" gauge links. That is, we integrate over all links on the boundaries of the 
lattice |l^. This is motivated by the fact that the gauge field for a uniform magnetic field 
in the usual continuum theory is not periodic [cf. A'f^^ = xHSn2, which generates Eq. (0)]. 
In a perturbative calculation moreover, the continuum external potential induces quantum 
fluctuations which are also aperiodic. Free boundary conditions seem most appropriate to 
describe the corresponding physics on the lattice. Free boundary conditions also serve to 
eliminate long-lived metastable states that occur in the U(l) theory with periodic boundary 
conditions, corresponding to closed Dirac strings winding through the (dual) lattice [|T3|. In a 
simulation of the U(l) external fleld problem on a lattice with periodic boundary conditions, 
it is necessary to use a global updating procedure to eliminate these strings ||13|- 

We compute the energy S{H) in the "induced" gauge flelds in lattice units according to 
(Refs. [|l|,|lU|,|ni) 
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and L is the length of a side of the lattice. To get the physical energy £^phys we perform 
a vacuum subtraction, and convert to units of the coupling constant g (which in three 
dimensions has units of (mass)^/^). For SU(2): 



^phys=/(i/9)'p(^)-^(0) 



(13) 



We now present results of our simulations of the U(l) and SU(2) theories. We have 
performed an extensive set of calculations over a wide range of values of H, f3, and lattice 
sizes. We flnd that relatively modest statistics (similar to those employed in the previous 
four-dimensional studies of Refs. II^;!!]) are sufficient to identify the main features of the 
effective potential. We use a 10-hit Metropolis algorithm. Our main results in the SU(2) case 
were obtained on a 32^ lattice, using 1000 sweeps to thermalize at each value of H and j3, 
followed by 5000 sweeps, keeping only every fourth conflguration for data (this results in an 



integrated autocorrelation time in the energy data which in most cases satisfies Tint ^ 0.6). 
To check finite volume effects, we have also done some calculations in the SU(2) case on a 
16^ lattice with an eightfold increase in statistics. 

Our U(l) calculations provide a useful check of our SU(2) code — we simply identify 
the link variables f/!J^^'^ = to do the U(l) simulation using our SU(2) code. Figure |1| 
shows results obtained on a 16^ lattice (using 1000 sweeps to thermalize at each value of H, 
followed by 10000 sweeps, keeping only every fourth configuration for data). The agreement 
with the expected classical behavior in the energy and in the induced field is very good, to 
within about 5% at /3 = 10 (the difference from the classical values is due mainly to finite 
size effects, and decreases with increasing (3). 

Figures ^ and || show the energy in the SU(2) theory on a 32^ lattice as a function of if, 
for various values of j3. We note that at large H one must run at comparatively large values 
of j3 in order to approach the scaling limit, as evident in Fig. ^ [Roughly speaking, we may 
expect some components of the "induced" field to be of 0{H), corresponding to a plaquette 
angle of 0{a?gH), which should remain ^ 1 in order to approach the continuum limit (cf. 
Ref. [|10|)-] Hence there may be appreciable finite volume effects at the largest fields in our 
data. 

In the region H/ g^ ^ 1 on the other hand, our data shows both good scaling behavior (cf. 
Fig. ^), and negligible finite volume effects, as shown in Fig. ^ where we compare our results 
from 16^ and 32^ lattices. Our data at /3 = 7 and 10 show clear evidence of a minimum in 
the energy. The data at /? = 7 lie about eight standard deviations below zero in a region of 
H/ g^ near 1. [The existence of a minimum is somewhat ambiguous in the case of our data 
at /3 = 12, since the errors in £^phys with the present level of statistics are large compared 
with its magnitude for /5 = 12 and H/ g^ < 2 (finite volume effects in these data are also 
significant, cf. Fig. ^). However, the data at /3 = 12 are generally consistent within errors 
with the energy obtained at /3 = 7 and 10.] The energy is generally more than an order of 
magnitude smaller than the classical result ^if^ over the whole range of H that we have 
studied. 

In summary, we have developed a method for the introduction of an external field in 
nonAbelian lattice gauge theories that follows from a coupling to an external current, as 
used in continuum field theories. We have performed a systematic study of scaling and 
finite volume effects, and we have obtained results for the effective potential over a wide 
range of field strengths. We have obtained the first evidence from lattice simulations of the 
existence of a nontrivial minimum in the effective potential for a background color-magnetic 
field in nonAbelian gauge theory. Our results support the Savvidy "ferromagnetic" picture of 
gluon condensation in the QCD vacuum. In future work, we will use the methods developed 
here to compute the effective potential in four-dimensional lattice SU(2). It should also 
be possible to look for the nonperturbative vortex structure associated with the "would- 
be" Nielsen- Olesen unstable modes ||^J§]. The three-dimensional gauge theory studied here 
should provide a very convenient system for such a study. 
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FIGURES 

FIG. 1. Effective potential and expectation value of the induced field strength for 
three-dimensional U(l) lattice gauge theory, as functions of the applied magnetic field strength 
H. These results are from a 16^ lattice at /3 = 10. 

FIG. 2. Effective potential for three-dimensional SU(2) lattice gauge theory, as a function of 
the applied color-magnetic field strength H, for various values of /?. These results were obtained 
on a 32^ lattice. 

FIG. 3. Effective potential for three-dimensional lattice SU(2), as in Fig. 0, here shown on an 
expanded scale for "small" field strengths H/g"^ < 2. 

FIG. 4. Effective potential for three-dimensional lattice SU(2) on two different lattice sizes. 
Data points for H/g^ < 1 were calculated at /3 = 7, while the points at larger H were taken at 
/3 = 12. 



